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Abstract 

In many important problems, one uses the median instead of the mean to estimate 
a population's center, since the former is more robust. But in general, computing the 
median is considerably slower than the standard mean calculation, and a fast median 
algorithm is of interest. The fastest existing algorithm is quickselect. We investigate a 
novel algorithm binmedian, which has 0(n) average complexity. The algorithm uses a 
recursive binning scheme and relies on the fact that the median and mean are always at 
most one standard deviation apart. We also propose a related median approximation 
algorithm binapprox, which has 0(n) worst-case complexity. These algorithms are highly 
competitive with quickselect when computing the median of a single data set, but are 
significantly faster in updating the median when more data is added. 

1 Introduction 

In many applications it is useful to use the median instead of the mean to measure the 
center of a population. Such applications include, but are in no way limited to: biological 
sciences; computational finance; image processing; and optimization problems. Essentially, 
one uses the median because it is not sensitive to outliers, but this robustness comes at 
a price: computing the median takes much longer than computing the mean. Sometimes 
one finds the median as a step in a larger iterative process (like in many optimization 
algorithms), and this step is the bottleneck. Still, in other problems (like those in biological 
applications), one finds the median of a data set and then data is added or subtracted. In 
order to find the new median, one must recompute the median in the usual way, because 
the standard median algorithm cannot utilize previous work to shortcut this computation. 
This brute-force method leaves much to be desired. 

We propose a new median algorithm binmedian, which uses buckets or "bins" to recur- 
sively winnow the set of the data points that could possibly be the median. Because these 
bins preserve some information about the original data set, the binmedian algorithm can use 
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previous computational work to effectively update the median, given more (or less) data. 
For situations in which an exact median is not required, binmedian gives rise to an approx- 
imate median algorithm binapprox, which rapidly computes the median to within a small 
margin of error. We compare binmedian and binapprox to quickselect, the fastest existing 
algorithm for computing the median. We first describe each algorithm in detail, beginning 
with quickselect. 

2 Quickselect 

2.1 The Quickselect Algorithm 

Quickselect (Floyd & Rivest 1975a) is not just used to compute the median, but is a more 
general algorithm to select the /c th smallest element out of an array of n elements. (When 
re is odd the median is the /c th smallest with k = (n + l)/2, and when re is even it is the 
mean of the elements k = n/2 and k = re/2 + 1.) The algorithm is derived from quicksort 
(Hoare 1961), the well-known sorting algorithm. First it chooses a "pivot" element p from 
the array, and rearranges ("partitions") the array so that all elements less than p are to 
its left, and all elements greater than p are to its right. (The elements equal to p go on 
either side). Then it recurses on the left or right subarray, depending on where the median 
element lies. In pseudo-code: 

Quickselect (^4, k): 

1. Given an array A, choose a pivot element p 

2. Partition A around p; let A\, A2, A3 be the subarrays of points <, =, > p 

3. If k < len(^4i): return Quickselect (A\ , k) 

4. Else if k > len(yli) + len(A 2 ): return Quickselect^, k - len(Ai) - len(^ 2 )) 

5. Else: return p 

The fastest implementations of quickselect perform the steps 3 and 4 iteratively (not 
recursively), and use a standard in-place partitioning method for step 2. On the other 
hand, there isn't a universally accepted strategy for choosing the pivot p in step 1. For the 
comparisons in sections [5] and El we use a very fast Fortran implementation of quickselect 
taken from Press, Teukolsky, Vetterling & Flannery (1992). It chooses the pivot p to be the 
median of the bottom, middle, and top elements in the array. 
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2.2 Strengths and Weaknesses of Quickselect 



There are many advantages to quickselect. Not only does it run very fast in practice, but 
when the pivot is chosen randomly from the array, quickselect has 0(n) average computa- 
tional complexity (Floyd & Rivest 19756). Another desirable property is that uses only 0(1) 
space, meaning that the amount of extra memory it needs doesn't depend on the input size. 
Perhaps the most underappreciated strength of quickselect is its simplicity. The algorithm's 
strategy is quite easy to understand, and the best implementations are only about 30 lines 
long. 

One disadvantage of quickselect is that it rearranges the array for its own computational 
convenience, necessitated by the use of in-place partitioning. This seemingly innocuous side- 
effect can be a problem when maintaining the array order is important (for example, when 
one wants to find the median of a column of a matrix). To resolve this problem, quickselect 
must make a copy of this array, and use this scratch copy to do its computations (in which 
case it uses 0(n) space). 

But the biggest weakness of quickselect is that it is not well-suited to problems where 
the median needs to be updated with the addition of more data. Though it is able to 
efficiently compute the median from an in-memory array, its strategy is not conducive to 
saving intermediate results. Therefore, if data is appended to this array, we have no choice 
with quickselect but to compute the median in the usual way, and hence repeat much of 
our previous work. This will be our focus in section [61 where we will give an extension of 
binmedian that is able to deal with this update problem effectively. 

3 Binmedian 

3.1 The Binmedian Algorithm 

First we give a useful lemma: 

Lemma. If X is a random variable having mean \i, variance a 2 , and median m, then 
m G [(j, — a, fj, + a] . 

Proof. Consider 



H — m\ = 



\E(X - m)\ 



< 



< 



E\X - m 
E\X-n\ 



< 
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Here the first inequality is due to Jensen's inequality, and the second inequality is due to 
the fact that the median minimizes the function ^(a) = ~E\X — a\. The third inequality is 
due to the concave version of Jensen's inequality. □ 

Given data points X\,...x n and assuming that n is odd3, the basic strategy of the 
bin median algorithm is: 

1. Compute the mean /x and standard deviation a 

2. Form B bins across \p — a, [i + a] , map each data point to a bin 

3. Find the bin b that contains the median 

4. Recurse on the set of points mapped to b 

Here B is a fixed constant, and in step 2 we consider B equally spaced bins across [/i — 
cr, [i + a]. That is, we consider intervals of the form 

' i i + l \ 

M - o- + -g • 2cr, — — • 2a J 

for i = 0, 1, . . . B — 1. Iterating over the data xi, . . . x n , we count how many points lie in 
each one of these bins, and how many points lie to the left of these bins. We denote these 
counts by Ni and Nl. 

In step 3, to find the median bin we simply start adding up the counts from the left 
until the total is > (n + l)/2. That is, b is minimal such that 

N L + £N,>"-±i. 

i=0 

Note that our lemma guarantees that b £ {0, . . . B — 1}. 

Finally, in step 4 we iterate through x\, . . . x n once again to determine which points were 
mapped to bin b. Then we consider B equally spaced bins across \i — a + [-jj • 2a, ■ 2a) , 
and continue as before (except that now Nl is the number of points to the left of the 
median bin, i.e. Nl <— Nl + Y^i=o ^i)- We stop when there is only one point left in the 
median bin - this point is the median. (In practice, it is actually faster to stop when the 
remaining number of points in the median bin is < C, another fixed constant, and then find 
the median directly using insertion sort.) 

An analogous strategy holds when n is even. Now at each stage we have to keep track of the bin that 
contains the "left" median (the k smallest element with k = n/2) and the bin that contains the "right" 
median (k — n/2 + 1). This is more tedious than the odd case, but is conceptually the same, so hereafter 
we'll assume that n is odd whenever we talk about binmedian. 
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Notice that a larger value of B gives fewer points that map to the median bin, but it 
also means more work to find b at each recursive step. For the time trials in sections [5] 
and[f3 we choose B = 1000 because it weighs these two factors nicely. We also replace the 
recursive step with an iterative loop, and use C = 20. Binmedian is implemented in Fortran 
for the time trials, and this code, as well as an equivalent implementation in C, is freely 



available at http://stat.stanford.edu/~ryantibs/median/ 



3.2 Strengths and Weaknesses of Binmedian 

Beginning with its weaknesses, the binmedian algorithm will run faster with fewer points 
in the median's bin, and therefore its runtime is necessarily dependent on the distribution 
of the input data. The worst case for binmedian is when there is a huge cluster of points 
surrounding the median, with a very large standard deviation. With enough outliers, the 
median's bin could be swamped with points for many iterations. 

However, we can show that the algorithm's asymptotic running time is largely unaffected 
by this problem: with very reasonable assumptions on the data's distribution, binmedian 
has 0(n) average complexity (see the Appendix). Moreover, this robustness also seems to 
be true in practice. As we shall see in section the algorithm's runtime is only moderately 
affected by unfavorable distributions, and in fact, binmedian is highly competitive with 
quickselect in terms of running time. 

Like quickselect, the binmedian algorithm uses O(l) space, but (like quickselect) this 
requires it to change the order of the input data. The greatest strength of binmedian is that 
the algorithm is able to quickly cache some information about the input data, when it maps 
the data points to bins. This gives rise to a fast median approximation algorithm, binapprox. 
It also enables binmedian, as well as binapprox, to efficiently recompute the median when we 
are given more data. We'll see this in section [6l but first we discuss binapprox in the next 
section. 



4 Binapprox 

4.1 The Binapprox Algorithm 

In some situations, we do not need to compute the median exactly, and care more about 
getting a quick approximation. Binapprox is a simple approximation algorithm derived from 
binmedian. It follows the same steps as binmedian, except that it stops once it has computed 
the median bin b, and just returns the midpoint of this interval. Hence to be perfectly clear, 
the algorithm's steps are: 

1. Compute the mean fi and standard deviation a 

2. Form B bins across \p — a, \i + a] , map each data point to a bin 
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3. Find the bin b that contains the median 

4. Return the midpoint of bin b 

The median can differ from this midpoint by at most 1/2 the width of the interval, or 
1/2 • 2a/B = <j/B. Since we use B = 1000 in practice, binapprox is accurate to within 
l/1000 th of a standard deviation. 

Again, we implement binapprox in Fortran for sections [5] and [6l This code as well as C 
code for binapprox is available at |http://stat.stanford.edu/~ryantibs/median/[ 

4.2 Strengths and Weaknesses of Binapprox 

Unlike binmedian, the runtime of binapprox doesn't depend on the data's distribution, since 
it doesn't perform the recursive step. It requires O(l) space, and doesn't perturb the input 
data, so rearranging is never a problem. The algorithm has 0(n) worst-case computational 
complexity, as it only needs 3 passes through the data. It consistently runs faster than 
quickselect and binmedian in practice. Most importantly it extends to a very fast algorithm 
to deal with the update problem. 

Binapprox's main weakness is fairly obvious: if the standard deviation is extremely large, 
the reported approximation could be significantly different from the actual median. 

5 Speed Comparisons on Single Data Sets 

We compare quickselect, binmedian, and binapprox across data sets coming from various 
distributions. For each data set we perform a sort and take the middle element, and report 
these runtimes as a reference point for the slowest way to compute the median. 

We test a total of eight distributions. The first four — uniform over [0, 1], normal with 
(j, = 0, a = 1, exponential with A = 1, chi-square with k = 5 — are included because they are 
fundamental, and frequently occur in practice. The last four are mixed distributions that 
are purposely unfavorable for the binmedian algorithm. They have a lot of points around 
the median, and a huge standard deviation. Two of these are mixtures of standard normal 
data and uniformly distributed data over [-10 3 ,10 3 ], resp. [-10 4 , 10 4 ]. The other two are 
mixtures of standard normal data and exponential data with A = 10 -3 , resp. A = 10 . 
All these mixtures are even, meaning that an equal number of points come from each of the 
two distributions. Table Q] gives the results of all these timings, which were performed on 
an Intel Xeon 64 bit, 3.00GHz processor. 

For the first four distributions, quickselect and binmedian are very competitive in terms of 
runtime, with quickselect doing better on the data from £7(0, 1) and N(0, 1), and binmedian 
doing better on the data from E(l) and x|- It makes sense that binmedian does well on the 
data drawn from E(l) and x|, upon examining the probability density functions of these 
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98.45 (1.13) 


1.59 
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Binmedian 


88.90 (1.63) 
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Table 1: Runtimes in milliseconds of quickselect, binmedian, binapprox, and sort for com- 
puting the median of 1 + 10 7 data points from different distributions. We performed 10 
repetitions for each distribution, and report the average and standard deviation of the run- 
times. (In fact, for each repetition we drew 20 different data sets from the distribution, and 
timed how long it takes for the algorithm to sequentially compute the medians of these 20 
data sets, in one contiguous block. Then we divided this number by 20; this was done to 
avoid timing inaccuracies.) We also list the runtime ratios, relative to the fastest time. 



7 



distributions. The density function of the exponential distribution, for example, is at its 
highest at and declines quickly, remaining in sharp decline around the median (log 2/ A). 
This implies that there will be (relatively) few points in the median's bin, so binmedian gets 
off to a good start. 

Generally speaking, binmedian runs slower than quickselect over the last four distribu- 
tions. But considering that these distributions are intended to mimic binmedian's degenerate 
case (and are extreme examples, at that), the results are not too dramatic at all. The bin- 
median algorithm does worse with the normal/uniform mixtures (skewed on both sides) 
than it does with the normal/exponential mixtures (skewed on one side). At its worst, 
binmedian is 1.3 times slower than quickselect (130.9 versus 101 ms, on the N(0, 1) and 
U(— 10 4 , 10 4 ) mixture), but is actually slightly faster than quickselect on the N(0, 1) and 
.E^IO -3 ) mixture. 

Binapprox is the fastest algorithm across every distribution. It maintains a reasonable 
margin on quickselect and binmedian (up to 1.62 and 1.69 times faster, respectively), though 
the differences are not striking. In the next section, we focus on the median update problem, 
where binmedian and binapprox display a definite advantage. 

6 Speed Comparisons on the Update Problem 

Suppose that we compute the median of some data x±, . . . x no , and then we're given more 
data x no +i, . . . x n , and we're asked for the median of the aggregate data set x±, . . . x n . Of 
course we could just go ahead and compute the median of xi, . . . x n directly, but then we 
would be redoing much of our previous work to find the median of xi, . . . x no . We'd like a 
better strategy for updating the median. 

Consider the following: we use binmedian to compute the median of x±, . . . x no , and save 
the bin counts Ni, % = 0, . . . B — 1, and Nl, from the first iteration. We also need save the 
mean [1q and the standard deviation cjq. Given x no +i, . . .x n , we map these points to the 
original B bins and just increment the appropriate counts. Then we compute the median 
bin in the usual way: start adding Nl + Nq + Ni + . . . and stop when this is > (n + l)/2. 

But now, since we haven't mapped the data X\, . . . x n to bins over its proper mean and 
standard deviation, we aren't guaranteed that the median lies in one of the bins. In the 
case that the median lies to the left of our bins (i.e. Nl > (n + l)/2) or to the right of our 
bins (i.e. Nl + Nq + . . . Nb-i < (n + l)/2), we have to start all over again and perform the 
usual binmedian algorithm on xi, . . . x n . But if the median does lie in one of the bins, then 
we can continue onto the second iteration as usual. This can provide a dramatic savings in 
time, because we didn't even have to touch xi, . . . x no in the first iteration. 

We can use binapprox in a similar way: save ctq, Ni, Nl from the median computation 
on xi, . . . x no , and then use them to map x no+ i, . . . x n . We then determine where the median 
lies: if it lies outside of the bins, we perform the usual binapprox on xi, . . . x n . Otherwise 
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we can just return midpoint of the median's bin. 

In this same situation, quickselect does not have a way to save its previous work, and must 
redo the whole computation^ In tabled we demonstrate the effectiveness of binmedian and 
binapprox on this update problem. Instead of adding data x no+ i, . . . x n only once, we add it 
many times, and we record how long it takes for the algorithms to compute these medians 
in one contiguous block. Binmedian and binapprox each keep a global copy of /io, o"0j -^i> Nl. 
They set these globals given the initial data x\,.. .x no , and as more data is added, they 
dynamically change them if the median ever lies outside of the bins. 

The first two situations from table [2] represent the best case scenario for binmedian and 
binapprox. The added data sets are fairly small when compared to the initial data set 
(n — no = 10 5 data points versus no = 1 + 10 data points). Also, the added data lies 
mainly inside the initial bin range [jUo — ^Ch "o + coL so the median will never lie outside of 
the bins. This means that binmedian and binapprox never have to undo any of their previous 
work (i.e. they never have to recompute the mean or standard deviation of the data set, 
and re-bin all the points), and are able update the median very quickly. At their best, 
binmedian is faster than quickselect by a factor of 3.97, and binapprox is faster by a factor 
of 24.55. 

The last two situations represent realistic, not-so-ideal cases for binmedian and binapprox. 
The added data sets are equal in size to the initial data set (n — no = 10 6 and no = 1 + 10 6 ). 
Moreover, the added data lies mainly outside of the initial bin range, which means that the 
median will eventually lie outside of the bins, forcing binmedian and binapprox to recompute 
the mean and standard deviation and re-bin all of the data points. (This may happen more 
than once.) Nonetheless, binmedian and binapprox still maintain a very good margin on 
quickselect, with binmedian running up to 3.94 times faster than quickselect and binapprox 
running up to 19.01 times faster. 

It is important to note that when binapprox successfully updates the median (that is, it 
finds the new median bin without having to recompute the mean, standard deviation, etc.), 
it is accurate to within cro/1000, where oq is the standard deviation of the original data set 
xi, . . . x ri0 . To bound its accuracy in terms of the standard deviation a of the whole data 

One might argue that quickselect sorts half of the and it could use this sorted order when 

given x„ +i, . . . x n to expedite the partitioning steps, and thus quickly compute the median of Xi,...x n . 
First, keeping this sorted order does not actually provide much of a speed up. Second, and most importantly, 
that quickselect might keep xi, . . . x no in the order that it produced is an unrealistic expectation in practice. 
There may be several other steps that occur between initially computing the median and recomputing the 
median, and these steps could manipulate the order of the data. Also, if the data's order shouldn't be 
perturbed in the first place (for example, if we're computing the median of a column of a matrix), then 
quickselect has to make a copy of x\, . . . x no , and only rearranges the order of this copy. It's not practical to 
keep this copy around until we're given x no +i, ■ ■ -x n , even if no is only reasonably large. 
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Binapprox 


125 (5.68) 


1 



Table 2: Runtimes in milliseconds of quickselect, binmedian, and binapprox for the median 
update problem. The algorithms initially compute the median of a data set, then we sequen- 
tially add more data, and the algorithms must recompute the median after each addition. 
Whenever new data is added, quickselect computes the median of the aggregate data set 
with its usual strategy, and hence repeats much of its previous work. On the other hand, 
binmedian and binapprox can use their cached bin counts to effectively update the median. 
For each distribution, we performed 10 repetitions, and we list the average and standard 
deviation of these times. We give the runtime ratios, which are taken relative to the fastest 
time. 
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set x\, . . . x n , consider 
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where is the mean of x±, . . . x, 
error in this situation is at most 



"0 



and [i is the mean of xi, . . . x n . Therefore, binapprox's 



0o 



< 



— a 

no 



(1) 



1000 - 1000 

If n = 2no, this becomes \/2<7/1000 ~ cr/700, which is still quite small. 

Finally, observe that we can use an analogous strategy for binmedian and binapprox 
to recompute the median after removing (instead of adding) a subset of data. Suppose 
that, having computed the median of x\, . . . x n , we wish to know median when x no+ i, . . . x n 
are removed from the data set. Then as before, we save fi,a,Ni,NL from our median 
computation on xi, . . . x n , map x m+ i, . . .x n to bins, and decrement (instead of increment) 
the counts appropriately. 



7 Discussion and Summary 

We have proposed an algorithm binmedian, which computes the median by repeated binning 
on a selectively smaller subset of the data, and a corresponding approximation algorithm 
binapprox. Binmedian has 0(n) average complexity when mild assumptions are made on the 
data's distribution function, and binapprox has 0(n) worst-case complexity. For finding the 
median of a single data set, binmedian is highly competitive with quickselect, although it 
becomes slower when there are many points near the median and the data has an extremely 
large standard deviation. On single data sets, binapprox is consistently faster no matter the 
distribution. Binmedian and binapprox both outperform quickselect on the median update 
problem, wherein new data is successively added to a base data set; on this problem, 
binapprox can be nearly 25 times faster than quickselect. 

In most real applications, an error of er/1000 (an upper bound for binapprox's error) 
is perfectly acceptable. In many biological applications this can be less than the error 
attributed to the machine that collects the data. Also, we emphasize that binapprox performs 
all of its computations without side effects, and so it does not perturb the input data in any 
way. Specifically, binapprox doesn't change the order of the input data, as do quickselect and 
binmedian. For applications in which the input's order needs to be preserved, quickselect 
and binmedian must make a copy of the input array, and this will only increase the margin 
by which binapprox is faster (in tables [T] and [2j) . 

Finally, it is important to note the strategy for binmedian and binapprox is well-suited 
to situations in which computations must be distributed. The computations to find the 
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mean and standard deviation can be easily distributed - hence with the data set divided 
into many parts, we can compute bin counts separately on each part, and then share these 
counts to find the median bin. For binapprox we need only to report the midpoint of this 
bin, and for binmedian we recurse and use this same distributed strategy. This could work 
well in the context of wireless sensor networks, which have recently become very popular in 
many different areas of research, including various military and environmental applications 
(see Zhao & Guibas (2004)). 

For any situation in which an error of cr/1000 is acceptable, there seems to be no 
downside to using the binapprox algorithm. It is in general a little bit faster than quickselect 
and it admits many nice properties, most notably its abilities to rapidly update the median 
and to be distributed. We hope to see diverse applications that benefit from this. 

Appendix 

Proof that Binmedian Has 0(n) Expected Running Time 

Let Xi, . . . X n be iid points from a continuous distribution with density /. Assuming that 
/ is bounded and that ~E(Xf) < oo, we show that the binmedian algorithm has expected 
0(n) running timejfl 

Now for notation: let / < M, and ~E(Xf) = a 2 . Let fi and a be the empirical mean 
and standard deviation of X±, . . . X n . Let Nj be the number of points at the j th iteration 
(Nq = n), and let T be the number of iterations before termination. We can assure T < n@ 

At the j th iteration, the algorithm performs Nj + B operations; therefore its runtime is 

T n n 

R = + B) < ^(Nj +B) = J2Nj + B(n + 1) 

j=0 j=0 j=0 

and its expected runtime is 

n 

E(R) < ^2'E(N j ) + B(n+l). 

3=0 

3 Notice that this expectation is taken over draws of the data X%, . . .X n from their distribution. This is 
different from the usual notion of expected running time, which refers to an expectation taken over random 
decisions made by the algorithm, holding the data fixed. 

4 To do this we make a slight modification to the algorithm, which does not change the running time. 
Right before the recursive step (step 4), we check if all of the data points lie in the median bin. If this is 
true, then instead of constructing B new bins across the bin's endpoints, we construct the bins across the 
minimum and maximum data points. This assures that at least one data point will be excluded at each 
iteration, so that T < n. 
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Now we compute E(iVj). Let b denote the median bin; this is a (random) interval of 
length 2a/ B. Well 

E(JVi) = nP(X 1 G 6) 

< nP(Xi G 6; a < a + e) + nP(a > a + e). 

Well P(Xi G b; a < a + e) < M ■ 2(a + e)/B (see the picture), and P(<r - a > e) = C/n 
(see the lemma; this is why we need a finite fourth moment of ^i). Therefore 

TL 

E(iVi) < -2M(a + e) + C. 
Jd 

In general 

E(Nj) = nP(X 1 eb,X 1 £b 2 ,... and X 1 G bj) = nP(X 1 G ^) 
where 6j is the median's bin at the i th iteration, so that we can similarly bound 

TL 

E(N J )<—2M(a + e)+C. 
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Therefore 



E(i?) < {w 2M ^ + e ) + °) + 5 ( n + ^ 

= 2M{o + e)\J2 V B * ) n+ (B + C)(n + 1) 



< 2M(a + e) ^ 1/B j j n + (£ + C){n + 1) 
\j'=o 

= O(n). 



Lemma. Wii/i i/ie same assumptions as above, P(<7 — cr > e) = 0(n 1 ) for any e > 0. 
Proof. Notice 

{a - a > e} = {(a - a) 2 > e 2 ; a > a) 

= {a 2 - 2(5" o + a 2 > e 2 ; a > a) 
C [a 2 - a 2 > e 2 ; <r > cr} 
= {<r 2 -a 2 >e 2 }. 

Recall that a 2 = ^ ^™ (Xj — X) 2 . Consider the unbiased estimator for a 2 , S 2 = ^rfO- 2 , and 

{d 3 -0 2 >€ i }c{lS 2 -0 2 >^}. 

By Chebyshev's inequality P(S* 2 — a 2 > e 2 ) < Var(5 2 )/e 2 . The rest is a tedious but 
straightforward calculation, using E(Xf) < oo, to show that Var(S' 2 ) = 0{n~ l ). □ 

Proof that Binmedian Has O(logn) Expected Number of Iterations 

Under the same set of assumptions on the data X\, . . . X n , we can show that the binmedian 
algorithm has an expected number of O(logn) iterations before it converges. We established 
that 



ft 

E(JVi) < - B 2M{a + e) + C, 



and by that same logic it follows 



Tfl 

~Ei(Nj\Nj-\ = m)< —2M(a + e) + C = m - g{m) 
B 
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where g{m) = m 2M{a+e / B 1] - C. By Theorem 1.3 of Motwani & Raghavan (1995), this 
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